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ABSTRACT 

As large optical quasar surveys for DLAs become a reality and the study of star forming gas in the 
early Universe achieves statistical robustness, it is now vital to identify and quantify the sources of 
systematic error. Because the nature of optically-selected quasar surveys makes them vulnerable to 
dust obscuration, we have undertaken a radio-selected quasar survey for DLAs to address this bias. 
We present the definition and results of this survey. We then combine our sample with the CORALS 
dataset to investigate the H I column density distribution function /hi(A'^, X) of DLAs toward radio- 
selected quasars. We find that fYii{N,X) is well fit by a power-law /Hi(iV, X) = kiN""^, with logfci 
= 22.90 and ai = — 2.18J^Q'2g. This power-law is in excellent agreement with that of optically-selected 
samples at low A^hi, an important yet expected result given that obscuration should have negligible 
effect at these gas columns. However, because of the relatively small size of the radio-selected sample, 
26 DLAs in 119 quasars, /hi(A^, X) is not well constrained at large iVni and the first moment of the H I 
distribution function, fig, is, strictly speaking, a lower limit. The power-law is steep enough, however, 
that extrapolating it to higher column densities implies only a modest, logarithmic increase in Qg. The 
radio-selected value of fig — l.lSlgggX 10"'^, agrees well with the results of optically-selected surveys. 
While our results indicate that dust obscuration is likely not a major issue for surveys of DLAs we 
estimate that a radio-selected sample of « 100 DLAs will be required to obtain the precision necessary 
to absolutely confirm an absence of dust bias. 

Subject headings: Galaxies: Evolution, Galaxies: 
Lines 
1. INTRODUCTION 



The introduction of large data sets from surveys such 
as the Sloan Digital Sky Survey (SDSS) has made possi- 
ble statistically significant studies of Damped Lya systems 
(DLAs ) (Prochaska, Herbert-Fort, & Wolfe 2005, here- 
after PHW05), quasar absorption systems defined as hav- 
ing an H I column density A^hi > 2 x 10^° cm~^ and which 
contain most of the neutral gas in the redshift interval 
z=[0, 5] (Wolfe, Gawiser, & Prochaska 2005). These large 
surveys for DLAs shed light on the history of the neutral 
gas content of the Universe and show how it is affected by 
star formation and gas replenishment. However, one ma- 
jor problem consistently affects magnitude-limited optical 
surveys: the issue of intervening dust and the possibility 
of obscuration bias. Because the metallicities of DLAs can 
be as high as 1/3 solar and are always above 1/1000 solar 
(Prochaska 2003), the presence of dust in these objects is 
not surprising. Specifically, evidence from element abun- 
dance patterns suggests the presence of depletion (Pettini 
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et al. 1999), while evidence for differential reddening sug- 
gests that dust obscuration is possible (Pei & Fall 1995). 
Dust obscuration in optically-selected surveys could be in- 
troducing a selection bias against DLAs whose high dust 
optical depth would hide the background quasar. This 
effect could seriously impact the results of statistics and 
derived values like Qg, the cosniological density of neutral 
gas (see Pei & Fall (1995)), particularly since high column 
density systems that dominate fi^, would be theoretically 
most likely to have strong dust obscuration. Here and 
throughout the paper we will use Qg to mean DP^^, the 
neutral gas contained in systems defined as being DLAs, 
i.e. with an H I surface density A'hi > 2 x 10^*^ cm~^. For 
a detailed discussion of the rationale behind this choice, 
see PHW05. 

Because radio observations are insensitive to the pres- 
ence of dust, a radio-selected sample of quasars does not 
suffer from this dust obscuration selection effect. There- 
fore, a radio-selected survey is a check on the possible in- 
troduction of biases in the magnitude-limited, optical sur- 
veys. One previous survey, the Complete Optical and Ra- 
dio Absorption Line System (CORALS) survey by Ellison 
et al. (2001), attempted to answer the question of the im- 
portance of dust obscuration by selecting quasars from the 
Parkes quarter- Jansky Flat-spectrum Radio Survey (Jack- 
son et al. 2002), and then following up with optical obser- 
vations of the selected quasars, regardless of magnitude, 
to search for the presence of DLAs. Ellison et al. found a 
slightly higher incidence of DLAs than that found by opti- 
cal surveys. From their measurement of Qg, Ellison et al. 
concluded that the effects of dust could be suppressing the 
magnitude-limited value of Qg by no more than a factor 
of two. 

The radio-selected surveys were motivated by several 
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studies which indicated that dust obscuration would signif- 
icantly bias the results of fig and other quantities derived 
from optical surveys. Fall & Pei (1993) constructed mod- 
els to predict the possibility of dust obscuration of quasars 
and found that between 10 and 70% of quasars at 2; = 3 
could be obscured, resulting in an underestimate of fig by 
the optical surveys. Recently, Wild & Hewett (2005) re- 
port on a survey of SDSS quasars for Ca II absorption-line 
systems with 0.84 < Zabs < 1-3. Using Ca II along with Fe 
and Mg lines, they claim that most of their sample likely 
contains DLAs and a significant amount of reddening. By 
modeling the reddening of these systems they make an es- 
timate that they are missing « 40% of Ca II systems from 
the SDSS due to dust obscuration, which they compare 
favorably to the upper limit of the CORALS survey. 

On the other hand, the SDSS-DR3 survey which found 
525 DLAs (PHW05), indicates that a dust bias, if present, 
is not an important effect. Prochaska et al. examine their 
results as a function of quasar magnitude, separating out 
the brightest 33% and the faintest 33% of the sample in 
each of four redshift bins. While there is not a statisti- 
cally significant difference in the line density, they mea- 
sure 40% higher ilg values towards brighter quasars at 
the 95% c.l. Since this is the opposite effect of what the 
dust bias would naively imply (a dust bias may imply 
that Q^'J'^sht ^ fj/'*™*^ since the brightest observed quasars 
should have less foreground dust obscuration which implies 
a smaller TVhi value and hence smaller fig), the SDSS-DR3 
results with a statistically significant higher value of fig to- 
wards brighter quasars point to the conclusion that dust 
obscuration is not an important effect. When Murphy 
& Liske (2004) examined the results of the SDSS-DR2, 
a sample including 70 DLAs they found no evidence for 
reddening. After examining the much larger SDSS-DR3 
Murphy et al. (2005) find evidence for reddening, but at a 
very low level - the implied dust to gas ratio is less than 
0.02, where dust to gas ratio is defined relative to that 
of the Milky Way (see equation 7 of Wolfe et al. (2003)). 
However, with all of these studies it is important to re- 
member that optical samples cannot measure dust bias in 
objects so heavily extincted that they are missing from 
the samples. As larger optical surveys for DLAs become 
feasible, due to surveys like the SDSS, and the statistical 
uncertainties become < 5%, potential causes of system- 
atic uncertainties, such as dust obscuration, must be fully 
understood. 

In this paper we will present the results of a radio- 
selected quasar survey that was undertaken by our group. 
This UCSD sample is approximately the same size as the 
CORALS sample and presents a comparable assessment of 
dust obscuration. We will analyze the combined results, 
attempt to assess the H I column density distribution func- 
tion, fiii{N, X), and show that our results for Jig do not 
differ in a statistically significant way from the results of 
optically-selected surveys, therefore suggesting that dust 
obscuration is most likely not a major problem affecting 
optically-selected quasar samples for DLAs. 

The organization of this paper is as follows: In § 2 
we describe the UCSD quasar sample, DLA identification 
method and analysis processes. In § 3 we review the stan- 
dard DLA statistical analysis methods. We discuss the 
results of the UCSD, CORALS and combined samples in 



§ 4, as well our process for dealing with the Empty Fields. 
And finally, in § 5 we will compare our results with the 
most recent optical surveys. 

Throughout the paper we will use the following cosmo- 
logical parameters (Bennett et al. 2003): 51a = 0.7,51m = 
0.3, Ho = 70km s-i Mpc^^. 

2. THE UCSD SAMPLE 

The UCSD sample consists of 68 objects selected from 
the 411 sources that comprise the complete Caltech-Jodrell 
Bank radio catalogs, including the Pearson-Readhead sam- 
ple (PR), the Caltech-Jodrell Bank VLBI Surveys 1 and 2 
(CJl and CJ2), and the Caltech-Jodrell Bank Flat-spectrum 
sample completion. While each sample was selected by a 
progressively lower flux density threshold, the sources for 
each sample were all selected to have declination (B1950) 
5 > 35° and Galactic latitude | 6 | > 10°. The PR sam- 
ple includes 65 objects with flux density at 6 cm (4859 
MHz), Seem > 1-3 Jy (Pearson & Readhead 1988). The 
CJl sample includes 135 sources with flux density at 6 
cm, 1.3Jy> Seem >0.7Jy (Polatidis et al. 1995; Thakkar 
et al. 1995; Xu et al. 1995), and the CJ2 sample consists 
of 193 mostly flat-spectrum objects with a flux density 
Seem >350 mJy (Taylor et al. 1994; Henstock et al. 1995). 
The CJF (CJF; Taylor et al. 1996) is a compilation of 
the flat-spectrum radio sources (spectral index flatter than 
o^uooMffz — —0.5) from the previous three surveys, plus 
and additional 18 sources for completion. 

An optical campaign to determine the type of source, 
magnitude, and redshift of the radio catalog objects fol- 
lowed, the results of which were compiled in the CJ cata- 
logs. The object optical identiflcation and determination 
of the R magnitude was done by automated scanning of 
the POSS plates or by eye. Redshifts were primarily taken 
from Vcron-Cctty, M.-P. & Veron, P. (1993) and from Hen- 
stock et al. (1997). 

From this large radio sample, our selection criterion in- 
cluded all objects identified as quasars with Zem > 2.0 
regardless of optical magnitude. The Zem > 2.0 cutoff was 
chosen to ensure sufficient spectral coverage to search for 
DLAs at wavelengths redward of the atmospheric cutoff 
at « 32OOA. We also included in our sample of 68 objects 
all 14 sources designated as optical empty fields in order 
to be sure that we were not artificially selecting objects 
brighter than an arbitrary optical magnitude. And finally, 
we included the 5 sources for which there was a tentative 
optical identification, but no redshift information. 

2.1. Observations and Analysis 

Observations of most quasar candidate objects were first 
carried out at Palomar, with follow-up done at Keck for 
faint or "Empty Field" (EF) objects (see Table 1). The 
majority of our spectra have better than 6 A FWHM. Five 
of the 68 objects in our sample were previously observed 
at moderate spectral resolution and the data or results 
for these objects were taken from the literature. These 
included quasars Q0014+813, Q0201+365, Q0636-t-680, 
Q0642-|~449 which were observed by Sargent et al. (1989), 
and quasar Ql 124+571 which was taken from Lanzetta et 
al. (1991). The iVni measurement for the DLA towards 
Q020H-365 was taken from Lu et al. (1993). 

Initial observations of most other targets were made 
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Table 1 
DETAILS OF OBSERVATIONS 



Ibicscopc 


Date 


Resolution 


No. of 




(No. nights) 


(at 4000 A) 


quasars observed 


Paloinar 


Nov 1995 (2) 


4-6 A 


17 


Palomar 


Apr 1996 (1) 


4-6 A 


8 


Palomar 


May 1996 (2) 


4-6 A 


9 


Palomar 


Doc 1996 (2) 


4-6 A 


16 


Paloinar 


June 1997 (1) 


4-6 A 


7 


Kock LRIS 


Nov 2001 (1) 


4-6 A 


3 


Keck ESI 


Apr 2002 (1) 


0.5 A 


14 


Keck ESI 


Aug 2002 (1) 


0.5 A 


7 


Keck LRIS 


Dec 2002 (1) 


11 A 


1 


Keck LRIS 


May 2003 (1) 


3-5 A 


14 



with the 200-inch Hale telescope of the Palomar Obser- 
vatory. Observations were made with the Double Spec- 
trograph and utilized gratings that were 300 lines mm~^ 
in the blue and the 315 lines mm^^ in the red, resulting 
in 4-6 A resolution using the l"slit, and wlO A resolution 
using the 2"slit when conditions were poor. All data were 
reduced using standard IRAF reduction packages. 

Follow up observations of empty fields and faint quasars 
were done at Keck with the Echellette Spectrograph and 
Imager (ESI, Sheinis et al. 2002) and the Low Resolu- 
tion Imaging Spectrometer (LRIS, Oke et al. 1995). With 
LRIS, the slit size was generally 0.7". The ESI observa- 
tions utilized the 1.0"slit in LowD mode and the 0.5"slit 
in echellette mode. 

Fifteen of the original 68 objects were excluded from the 
final statistical data set for the following reasons: Seven 
were determined to be galaxies, stars, or low redshift quasars, 
and the remaining eight were deemed empty fields (EFs). 
Table 2 contains the details of the final 53 quasars used 
in the UCSD statistics, while Table 3 and Table 4 contain 
a summary of the empty fields and the discarded objects 
respectively. 



Table 4 
DISCARDED OBJECTS 



Object 


Reason lor Discard 


Survey 


0843+575 


galaxy 


CJ2 


1125+596 


quasar at z^m — 1-78 


CJ2 


1308+471 


galaxy 


CJ2 


1436+763 


star 


CJ2 


1809+568 


No significant emission feature 


CJ2 


2238+410 


Spectrum dubious 


CJ2 


2319+444 


quasar at z^m — 1-24 


CJ2 



2.2. Damped Lya systems 

The Palomar data were reduced using standard IRAF 
packages, while the Keck data were reduced using IDL re- 
duction software'*. The reduced quasar spectra were con- 
tinuum fitted and normalized and then analyzed to find 
regions in which the restframe equivalent width of an ab- 
sorption feature was > 5 A and located in a region of good 
signal to noise. The equivalent width of the spectrum was 
calculated and an equivalent width array was then ana- 
lyzed for regions that were greater than the 5 A restframe 

*http: //www. ucolick.org/~xavier/IDL 



cutoff, as explained by Wolfe et al. (1995). We searched 
all regions of the spectrum blueward of Lya emission, be- 
ginning with the lowest wavelength at which the error was 
below the restframe equivalent width threshold of > 5 A at 
the 5cr level. All candidate detections were then inspected 
by eye to determine if they were indeed DLAs. False de- 
tections were usually quite obvious to exclude as blended 
fines, Ly/3, etc. 

Nine DLAs were found, two of which, towards quasars 
Q0824-I-355 and Q1014-h615, were within 3,000 km s'^ of 
the Lya emission peak, and therefore considered "associ- 
ated" . Following the standard practice, we exclude these 
"associated" DLAs from the sample in order to insure that 
we are not detecting objects that are physically associated 
with the quasar. Discarding these two leaves a final seven 
DLAs to be included in the UCSD sample. 

2.2.1. Lya Fits 

The DLA systems were fitted with Voigt profiles using 
the IDL tooP X-fitdla which allows the user to interactively 
modify the Voigt profile and continuum placement. In all 
but one case, that of Q1239+376, the DLA redshift was 
constrained by the corresponding metal lines with errors as 
given in Table 2. In the case of Q1239+376, the metal lines 
were too weak for use in constraining the DLA redshift 
and we instead determined the best fit interactively by eye 
using x_fitdla. This method results in larger uncertainties 
for Zabs and TVhi- 




3450 3500 3550 3600 3650 3700 3750 3800 3850 



Fig. 1. — Voigt profile DLA fits for the six new DLA systems in 
the UCSD sample. The best Voigt profile is indicated by the green 
curve, surrounded by profiles with A'^hi displacements ± O.lSdex (or 
± 0.25 dex) in red. The continuum placement can be seen as the 
purple dotted line, while the error array is represented by the blue 
dotted line. 



For most of our sample DLAs a conservative estimate 
of the uncertainty in A^hi is 0.15 dex. However in one case, 

^http://www. ucolick.org/~xavier/IDL 
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Table 2 
UCSD SURVEY SAMPLE 



Quasar 


^em 


R mag 


Nm X 10^" cm-" 


^abs 


^min 


^max 


6cm Flux (Jy) 


Survey^ 


Q0014+813 


3.384 


15.9 






1.591 


3.340 


0.551 


CJ2 


Q0153+744 


2.338 


16.0 






1.568 


2.305 


1.794 


PR 


Q0201+365 


2.912 


17.5 


2.4^ 


2.461 


1.632 


2.873 


0.349 


CJ2 


Q0212+735 


2.367 


19.0 






1.742 


2.333 


2.444 


PR 


Q0604+728 


3.53 


20.3 






2.651 


3.485 


0.654 


CJ2 


Q0609+607 


2.710 


19.0 






1.650 


2.673 


1.059 


CJ2 


Q0620+389 


3.470 


20.0 






1.842 


3.425 


0.87 


CJl 


Q0627+532 


2.200 


18.5 






1.619 


2.168 


0.485 


CJ2 


Q0636+680 


3.174 


16.2 






1.591 


3.132 


0.499 


CJ2 


Q0642+449 


3.406 


18.5 






1.591 


3.362 


0.78 


CJl 


Q0727+409 


2.501 


17.0 






1.578 


2.466 


0.468 


CJ2 


Q0738+491 


2.318 


21.0 






1.928 


2.285 


0.352 


CJ2 


Q0749+426 


3.590 


18.1 






2.118 


3.544 


0.461 


CJ2 


Q0800+618 


3.04 


19.6 


3.16 ±0.15 


2.9603 ±0.0017 


2.234 


3.000 


0.981 


CJ2 


Q0803+452 


2.102 


19.6 






1.594 


2.071 


0.414 


CJ2 


Q0824+355 


2.249 


19.7 


2.0'^ 


2.2433 


1.655 


2.217 


0.746 


CJ2 


Q0833+585 


2.101 


18.0 






1.546 


2.070 


1.11 


CJl 


Q0836+710 


2.180 


16.5 






1.507 


2.148 


2.423 


PR 


Q0902+490 


2.690 


17.2 






1.550 


2.653 


0.547 


CJ2 


Q0917+449 


2.180 


19.0 






1.578 


2.148 


0.80 


CJl 


Q0930+493 


2.590 


18.4 






1.666 


2.554 


0.574 


CJ2 


Q1014+615 


2.800 


18.1 


2.5'^ 


2.7670 


2.263 


2.757 


0.631 


CJ2 


Q1053+704 


2.492 


18.5 






1.801 


2.457 


0.71 


CJl 


Q1124+571 


2.890 


18.0 






1.796 


2.851 


0.597 


CJ2 


Q1144+542 


2.201 


20.5 






1.632 


2.169 


0.88 


CJl 


Q1155+486 


2.028 


19.9 






1.632 


1.998 


0.445 


CJ2 


Q1214+588 


2.547 


19.5 






1.632 


2.512 


0.307 


CJ2 


Q1239+376 


3.818 


19.5 


2.0 ±0.15 


3.4082^= 


2.344 


3.770 


0.446 


CJ2 


Q1325+436 


2.073 


18.5 






1.549 


2.042 


0.533 


CJ2 


Q1333+459 


2.449 


18.5 






1.612 


2.414 


0.76 


CJl 


Q1337+637 


2.558 


18.5 






1.550 


2.522 


0.431 


CJ2 


Q1413+373 


2.360 


17.3 






1.607 


2.326 


0.383 


CJ2 


Q1421+482 


2.220 


18.9 






1.549 


2.188 


0.536 


CJ2 


Q1427+543 


2.980 


20.7 






2.331 


2.940 


0.718 


CJ2 


Q1435+638 


2.068 


15.0 






1.591 


2.037 


1.24 


CJl 


Q1526+670 


3.020 


17.1 






1.977 


2.980 


0.417 


CJ2 


Q1547+507 


2.169 


18.5 






1.582 


2.137 


0.74 


CJl 


Q1602+576 


2.858 


16.8 






1.630 


2.819 


0.351 


CJ2 


Q1624+416 


2.550 


22.0 






1.732 


2.515 


1.632 


PR 


Q1645+635 


2.380 


19.4 


3.55 ±0.15 


2.1253 ±0.0003 


1.536 


2.346 


0.444 


CJ2 


Q1745+624 


3.886 


18.3 






3.085 


3.837 


0.580 


CJ2 


Q1755+578 


2.110 


18.6 


25.1 ±0.15 


1.9698 ±0.0009 


1.630 


2.079 


0.455 


CJ2 


Q1758+388 


2.092 


17.8 






1.512 


2.061 


0.92 


CJl 


Q1834+612 


2.274 


17.6 






1.599 


2.241 


0.590 


CJ2 


Q1839+389 


3.094 


19.5 


5.0 ±0.15 


2.7746 ±0.0009 


1.911 


3.053 


0.476 


CJ2 


Q1850+402 


2.120 


17.9 


20.0 ±0.25 


1.9888 ±0.0058 


1.669 


2.089 


0.535 


CJ2 


Q2015+657 


2.845 


19.1 






2.734 


2.807 


0.500 


CJ2 


Q2017+745 


2.187 


17.9 






1.602 


2.155 


0.500 


CJ2 


Q2136+824 


2.350 


18.9 






2.002 


2.317 


0.509 


CJ2 


Q2255+416 


2.150 


20.9 






2.119 


2.119 


0.99 


CJl 


Q2259+371 


2.228 


20.4 






1.632 


2.196 


0.406 


CJ2 


Q2310+385 


2.181 


17.5 






1.630 


2.149 


0.484 


CJ2 


Q2356+385 


2.704 


18.6 






1.771 


2.666 


0.449 


CJ2 



"^HI value taken from Lu ct al. (1993). 

''Associated systems: z^ts within ai 3,000 km s^-"- of Zem- 

'^Weak metals, therefore best fit determined by eye. 

'^CJl, CJ2 = Caltech Jodrell Bank, PR = Pearson- Readhead 
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Table 3 
EMPTY OR EXTENDED FIELDS 



Object RA(J2000) Dcc(J2000) Exp. Rh™' 

time (s) 3(7 



Morphology 



R mag 



6cm Flux 

(Jy) 



Survey 



0102+480 


01 05 49.93 


+48 19 03.19 


294 


26.1 


no detection 




1.080 


CJl 


0633+596 


06 38 02.87 


+59 33 22.21 


500 


26.4 


possibly extended 


25.8 ±0.7 


0.482 


CJ2 


0718+793 


07 26 11.74 


+79 11 31.0 


500 


26.1 


extended, «7"x3" 


24.5 ±0.4 


0.467= 


CJ2 


1107+607 


11 10 13.09 


+60 28 42.57 


600 


26.5 


extended, «4"x2.5" 


25.8 ±0.8 


0.400 


CJ2 


1205+544 


12 08 27.50 


+54 13 19.53 


600 


26.5 


possibly extended 


24.6 ±0.3 


0.397 


CJ2 


1312+533 


13 14 43.83 


+53 06 27.73 


600 


26.5 


extended, «3"x3" 


25.4 ±0.6 


0.433 


CJ2 


1828+399 


18 29 56.52 


+39 57 34.69 


900 


26.9 


no detection 




0.353 


CJ2 


2054+611'' 


20 55 38.84 


+61 22 00.64 


900 


26.7 


possibly extended 




0.414 


CJ2 



" Limiting magnitude per seeing element above sky background 

''R mag estimation of extended smudge 

'=@1.4GHz 

■^Uncertain identification, eitfier z = 1.588, 3.0, 3.3 



that of Q1850+402, where the Voigt profile proved diflicult 
to fit, we report an uncertainty of 0.25 dex in iVni- Figure 1 
shows Voigt profile fits for each DLA, except for the Zats = 
2.461 DLA towards Q0201+365, which is in the existing 
literature (Sargent et al. 1989). Now we will give brief 
details on each DLA system. 

Q0201+365: Nm fit taken from Lu et al. (1993) and 
discussion therein. 

Q0800+618: The difficulty in estimating the contin- 
uum placement in such close proximity to the Lya emis- 
sion peak made this DLA system somewhat diSicult to fit. 
There is also a possibility of some blended absorption. 

Q1239+376: High signal to noise and good placement. 
The only problem with this fit was some blending on the 
red side. 

Q1645+635: High signal to noise and lack of blending 
resulted in a good Voigt profile fit to this DLA profile. 

Q1755+578: Close proximity to Lya emission peak 
and blending on the blueward side made the fit difficult. 

Q1839+389: Straightforward fit and good continuum. 

Q1850+402: Close proximity to Lya emission peak, 
lower signal to noise and blending made this a more dif- 
ficult fit, yielding an increased error margin on the Voigt 
profile of 0.25 dex. 

3. DLA STATISTICS 

Our goal of determining the impact of dust obscuration 
in surveys of DLAs requires that we be able to compare 
our radio-selected survey to the results of optically-selected 
surveys. We will now introduce some of the standard sta- 
tistical quantities used to describe and quantify surveys of 
DLAs. 

3.1. Az, g(z) and n{z) 

The redshift path, Az, is defined as the total rcdshift 
interval along which a damped Lya feature with rest frame 
equivalent width exceeding 5A would be detected at more 
than 5(7 significance. It is defined as follows. 



n 

Az = ^{z^"'' - z™" 



(1) 



where the summation is over the n quasars in the survey, 
Zmin is determined to be the lowest spectral wavelength 
with good signal-to-noise, and z^ax is the redshift corre- 
sponding to the maximum spectral wavelength included in 
the search. We define Zmax by. 



= ZgsO - (1 + Zgso)/WO 



(2) 

This corresponds to 3000 km s~^ blueward of the Lya 
emission feature. This cutoff ensures that a damped Lya 
system is not physically associated with the quasar. 

The redshift path density, g{z), gives an idea of the sta- 
tistical significance as a function of redshift of a survey for 
DLAs. It is defined as the number of quasars with sight- 
lines containing a particular redshift interval (Lanzetta et 
al. 1991). Specifically, 



9 



iz) = j2H{zr''-z)H{z-zrn 



(3) 



where H is the Heaviside step function, the sum is over n 
quasars (Lanzetta et al. 1991) and. 



(4) 



Az — / g{z)dz 



where the integral is over all z paths in the survey. The 
DLA number density, n^z)^ is simply the number of DLAs 
per unit redshift. 



n{z) 



m 
A^ 
where m is the number of DLAs. 



(5) 



3.2. /hi(-/V, X) ; The H I Frequency Distribution 
Function 

Following the direction of previous works such as Lanzetta 
et al. (1991) and PHW05, we can define a neutral hydrogen 



The UCSD Radio- Selected Quasar Survey for DLAs 



frequency distribution function that describes the number 
of DLAs in a range of column densities, {N, N + dN), and 
a range of absorption distances, {X, X + dX), 

fmiN,X)dNdX , (6) 

where the absorption distance, AX, is defined as follows: 

Ho 



AX ^ dX = 



H{z) 



{l + zydz 



(7) 



where Hq is Hubble's constant. 



3.3. ^dla(^) •' The Damped Lya Line Density 

The zeroth moment of the H I frequency distribution 
function is known as the line-density of DLAs £dla(X). 
The line-density represents the number of systems per unit 
absorption distance and is defined as: 



tDLA 



(X) 



fmiN,X)dN 



(8) 



Nt 



As described in PHW05, the line-density is related to the 
covering fraction of DLAs on the sky. This relationship 
is apparent if we describe the frequency distribution func- 
tion in terms of an average cross-section A{X), and the 
comoving number density of DLAs nr)LAiX): 

Jm{N,X) = {c/Ho)nDLA{X)A{X) . (9) 

(see Wolfe, Gawiser, & Prochaska 2005, for details). 

3.4. rjg.- The Cosmological Neutral Gas Mass Density 

An important parameter in describing any quasar sur- 
vey for DLAs is the first moment of the H I frequency 
distribution function, the neutral gas mass density, Qg. 
It is believed that fig is closely related to the amount of 
neutral hydrogen available for star formation and hence, 
places an important tracer on the history of star forma- 
tion in the Universe. Acquiring this parameter through 
surveys for DLAs is an important constraint on the neu- 
tral gas reservoir available for star formation in the early 
(z > 2) Universe, fig is defined as follows: 



n^ix) EE 



^J.mHHo 



N„ 



CPc 



Nfm{N,X)dN (10) 



N„ 



where fj, is the mean molecular mass of the gas (taken to be 
1.3), m.H is the mass of the hydrogen atom, pc is the crit- 
ical mass density, fiii{N,X) is the frequency distribution 
function of neutral gas found in DLAs and the integration 
is from Nmin = 2 x 10^" cm^^ to N^ax = oo. We follow 
previous works, i.e. Lanzetta et al. (1991), and replace 
this frequency distribution function by its evaluation in 
the discrete limit as follows. 



ng = 



pmnHo I]iV(HI) 



CPc 



AX 



(11) 



where the sum is performed over the A'hi measurements 
of the DLA systems in a given redshift interval with total 
pathlength AX. As emphasized by PHW05 and discussed 



below, equation 11 only provides an accurate evaluation of 
equation 10 if the survey is sufficiently large that the ob- 
served fYn{N,X) distribution becomes steeper than 7V~^ 
at large A^hi- If this is not the case, equation 11 provides 
only a lower limit to Qg. 

4. RESULTS 

We will now describe the results of the UCSD radio- 
selected survey, the CORALS radio-selected survey, and 
the combination of these two surveys, which we will refer 
to as the combined sample. Details of the results of each 
survey are listed in Table 5. 

4.1. UCSD Survey Results 

The UCSD sample consists of 7 DLAs in 53 quasars of 
Zem > 2.0 with a total redshift path of Az = 41.15. This 
resulted in a number of DLAs per unit redshift, n{z) = 
0.171q'o7, where the error bars are the standard la Poisso- 
nian errors using Gehrels' tables for small number statis- 
tics (Gehrels 1986). Figure 2 presents g{z) versus z for 
the UCSD sample in green. The line density of DLAs over 
the cosmological redshift path of AX = 130.43 resulted 
in £dla(^) = 0.051q'q2i while the mass density of neutral 



gas is D,g = 0.84_ 



0.45 



X 10- 



While we report fig as a 



detection, it is strictly speaking, a lower limit because we 
do not measure fiii{N,X) to fall off faster than N~'^ (see 
§ 4.3). 
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Fig. 2. — Redshift sensitivity function g{z) as a function of red- 
shift for the UCSD survey, the CORALS survey, and the combined 
sample. 



4.2. CORALS Survey Results 

The Complete Optical and Radio Absorption Line Sur- 
vey (CORALS) (Elhson et al. 2001) was the first attempt 
to utilize a radio-selected quasar survey as a basis for a 
search for DLAs. They selected quasars from the complete 
Parkes quarter- Jansky flat spectrum sample (Jackson et 
al. 2002), comprised of 878 radio sources with spectral in- 
dex a2 7GHz > —0.4 and declinations between -1-2.5° and 
-80°. Ellison et al. (2001) hmitcd their data set to 66 
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Zem > 2.2 quasars in which they found 22 DLAs. Three 
of these DLAs were classified as "associated" and dropped 
from the final sample. Two more DLAs were excluded be- 
cause they fell outside of the range 1.8 < Zats < 3.5, that 
Ellison et al. (2001) set for their statistical sample, citing 
that there appears to be little evolution of f2g in this range. 
Since the UCSD sample did not have this ^^""^ cutoff, we 
included all 19 CORALS DLAs in the statistics for the 
combined sample. 

Detailed results of the CORALS survey, including the 
two DLAs in sightlines to quasars with Zabs > 3.5, are 
listed in Table 5. Statistics for the number density and 
neutral gas mass density in the CORALS survey resulted 
in their conclusion that previous, magnitude limited sur- 
veys could have underestimated these values by as much as 
a factor of two. The plot of g(z) versus z for the CORALS 
survey is shown in blue in Figure 2. Over a total redshift 
interval, Az = 57.16, the number of DLAs per unit red- 
shift, n{z) — G.33lg'|!)7. Over a cosmological redshift path 
of AX = 186.68, the Hue density of DLAs in the CORALS 
survey, £dla(X) = G.102j;[5:oi which is double that of the 
UCSD survey. CORALS compared their neutral gas mass 
density in DLAs Qg = l.STtoss ^ 10"^ '^ith the compila- 
tion of Peroux et al. (2001) and Rao & Turnshek (2000), 
and concluded that fig derived from optically-selected sur- 
veys could be underestimated by up to a factor of two. 

However, Ellison et al. do concede the uncertainty of 
their conclusion primarily because the small survey fails to 
fully sample the column density distribution and secondly 
because their high value of Qg is dominated by two rela- 
tively high column density systems (both incidentally in 
front of "moderately bright" quasars, B — 19.5, 20, which 
qualitatively matches the result of PHW05 that there is 
an anti-correlation between quasar magnitude and Nm). 

Table 5 
RESULTS 
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4.3. Combined Results 

For simplicity, we present detailed analyses for just the 
combined sample, which has the greatest statistical signif- 
icance. Figure 2 presents g{z) versus z for the combined 
sample in red. The CORALS sample begins abruptly at 
z = 1.8, the Zmin cutoff of their sample. The UCSD sam- 
ple continues down to a Zmin ~ 1-51 for some quasars. 
For z w 2 the combined sample is nearly double that 
of CORALS. For higher redshift intervals (i.e. z — 3), 
the CORALS survey contributes roughly 2/3 of the path- 
length. Of course, the combined sample gives the best con- 
strained estimate of the number density of n(z)= 0.261o;q5. 



The combined sample is large enough to attempt an 
analysis of the H I distribution function, /hi(A'^, X). This 
sample spans the redshift interval z — [1.51,4.4] with an 
integrated absorption pathlength /S.X ~ 317.11 and a col- 
umn density weighted mean redshift of 2.28. In practice, 
we can evaluate fYii{N,X) in the discrete limit and plot 
the resulting /hi(A^, X) in A'hi bins of some A7V. In Fig- 
ure 3, we plot in red /hi(A'^, X) for the combined sample, 
in A'hi bins of AiV = 0.4 dex, calculated in the following 
way: 



/Hl(iV,X) = 



AX 



(12) 



where ttidla is the number of damped Lya systems within 
{N, N + AN) in the AX interval and the error bars are de- 
termined by Poisson uncertainty at the 84% c.l. according 
to the value of moLA- Also plotted, in black, are the re- 
sults of the optically-selected SDSS-DR3 (PHW05) survey 
for comparison. Following PIIW05, we have ovcrplotted 
the best-fit solutions of two possible functional forms of 
/hi(A^, -^)- Because of the small sample size of this sur- 
vey we will attempt to fit only a single power-law and a 
F-function. The single power-law form is as follows: 



/Hi(iV,X)-fci7V"i , 
and the F-function is (e.g. Fall & Pei 1993): 



fm{N,X) 



N 



exp 



-N 



(13) 



(14) 



We have performed a maximum likelihood analysis to con- 
strain the parameters and set the constants fci and ^2- 
A summary of the fit parameters, along with those of 
the optically-selected SDSS-DR3 survey for comparison, 
is given in Table 6. The best fit slope of the single power- 



law is ai — —2.18 



+0.20 
-0.26- 



This single power-law slope can 



be compared favorably with the optical SDSS-DR3 survey 



single power-law slope of ai = — 2.19;too5 '-'^'^'" their en- 
tire redshift range, z = [2.2, 5.5]. This correspondence is 
expected because the radio-selected survey is dominated 
by the low column density end which matches that of the 
optical, and can be seen as a confirmation of the two tech- 
niques. 

While the single-power law gives a good fit to the radio- 
selected data, we also attempt to fit the F-function for 
the following two reasons: First, the single-power law is 
unphysical, i.e. the fit must turn over in order for Q,g to 
converge, and second, unlike the single-power law, the F- 
function provided a satisfactory fit to the optically-selected 
data. However, unlike the optically-selected sample, the 
radio-selected sample gives nearly the same fit for the F- 
function as for the single power-law, a2 = —2.Vlt_^'^. 
While we derive a formal value of the break column den- 



sity N..y — 22.49_Q'3g, we interpret this as an unrealistic 
extrapolation of the data. Rather, the small size of the 
radio-selected sample cannot reliably determine a break 
column density, and therefore we cannot demonstrate that 
our rjg converges. 

To determine if the radio-selected data rules out the 
optically-selected F-function fit, we performed a chi-squared 
test on the radio-selected data and optically-selected F- 
function fit. The results of the chi-squared test, Probi(x^ 
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Fig. 3.— The H I frequency distribution /hi(A^, ^) for the 26 
DLAs of the combined sample is plotted in red. Overplotted are 
the fits of a single power-law, the dot-dashed line in blue, and a F- 
function, the dashed line in red. The last bin contains the la upper 
limit. Plotted in black is the /hi(^i^) foi" the optical data from 
the SDSS-DR3, with the F-function fit in green. 



> 5.96) = 1.5% show that we can reject the fit at the 5% 
level, but not at the 1% level of confidence. While this may 
be evidence for modest disagreement between the two sam- 
ples, we interpret this disagreement to be primarily due to 
the fact that we cannot constrain the radio-selected fit at 
large iVni , i.e. the radio-selected sample does not contain 
enough DLAs to fully sample the H I distribution func- 
tion. We note that a more conservative two-sided KS test 
shows agreement between the radio-selected data and the 
optically-selected F-function fit at the 77% level. 



Table 6 
FITS TO /hi(A^,X) 



Form 
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Combined Radio^ 


Single 
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"^Prochaska, Herbert-Fort, .! 
''Mean absorption redshift : 
'^Mean absorption redshift : 



Wolfe (2005) 
: 3.06 
: 2.28 



The line density of DLAs in the combined sample, taken 
over the entire redshift interval, z = [1.5, 4.4], is ^dlaI-'^) = 
0.081q'q2 a-t a- median z = 2.35, where the errors represent 
the la Poisson uncertainty in ixidla- In Figure 4 we plot 
^dla(^) for the combined sample in red, evaluated in the 
discrete limit, 



tDLA 



(X) 



AX 



(15) 



We have grouped the data into four redshift bins, z — [1.5, 



2.2], [2.2, 2.5], [2.5, 3.0], and [3.0, 4.4] to allow for compar- 
ison with the results of PHW05. We have overplotted the 
results of the SDSS-DR3 optical survey in black. These 
points are grouped into redshift bins z = [2.2, 2.5], [2.5, 
3.0], [3.0, 3.5], [3.5, 4.0], [4.0, 5.3]. The black point marked 
by a star in redshift bin z — [1.5, 2.2] is a compilation of 
optical surveys for DLAs produced by Peroux et al. (2003). 
Although the central values of the line densities of the ra- 
dio and optically-selected surveys are different, the differ- 
ence is not statistically significant. Note the radio sample 
gives a somewhat higher line density at all redshifts, and it 
is interesting to note that beginning at z = 2.2, the trend 
of increasing line density with increasing redshift, z > 2.2, 
is present in both samples. In fact, the central values of 
the radio sample follow the same qualitative shape, even 
the unusual 'dip' at z « 2.3. While it would make sense 
that the SDSS-DR3 survey, with its statistically signifi- 
cant numbers of quasars and DLAs is actually detecting 
a physically meaningful trend - PHW05 claim the decline 
in ^dla(-^) is due to a decrease in DLA cross-section with 
time - the correspondence with the combined radio sam- 
ple, of relatively so few objects, is likely a coincidence. 
However, although the error bars are large, we can inter- 
pret this similarity in line density evolution with the sta- 
tistically significant results of the SDSS-DR3 as support 
of our results. 
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Fig. 4. — Plot of the line density of DLA systems £dla(^) versus 
redshift for the combined sample (solid lines). Overplotted is the 
^dla(^) for the optical data from the SDSS-DR3 survey and the 
Peroux compilation (dashed lines). 



While we report a detection of the neutral gas mass den- 



sity, rig = 1.15 



-1-0.37 
'0.38 



X 10 '^, our result is actually a lower 



limit due to the insufficient size of the combined sample. In 
Figure 5 we plot Ug for each of the UCSD, CORALS and 
combined samples. Errors are calculated using a modified 
bootstrap method as described by PHW05, and the val- 
ues are plotted at the iVni weighted mean redshift. Also 
plotted, in black, arc the Qg values determined by the 
optically-selected SDSS-DR3 survey (PHW05), which cov- 
ers a redshift range z > 2.2. Because the SDSS-DR3 sam- 
ple consists of over 500 DLAs the values of fig are plotted 
in five redshift bins: z = [2.2,2.5], [2.5, 3.0], [3.0, 3.5], [3.5, 
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4.0], and [4.0, 5.5]. And finally, plotted in the bin range 
z — [1.7,2.2], is the compilation by Peroux et al. (2003). 
It is seen that the lower limits of all of the radio-selected 
samples agree well with the optically-selected data. 
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Fig. 5. — The neutral gas mass density, 
CORALS and combined samples are plotted in green, blue and red 
respectively. For clarity, only the redshift bin of the combined sam- 
ple is plotted, Zcombined = [1-51, 4.41] (dotted red). The redshift 
bins of the radio samples (not plotted) are z^^^^ = [1.51, 3.84], and 
Zcorals = [1.80, 4.41] and the points are plotted at the JVhi weighted 
mean redshift. Also plotted are the values of Og as a function of red- 
shift for the optical SDSS-DR3 survey (black points). The Peroux 
compilation data point at z < 2.2 (marked with a cross) does not 
include measurements from the SDSS survey. All error bars are 1 

(T. 



4.4. Empty Fields 

The eight fields for which no optical identification of a 
quasar was obtained are called the "Empty Fields" (EFs). 
Table 3 contains the details of each EF while Figure 6 
contains Keck images of the EFs. These eight fields were 
determined to contain either nothing of significance, or 
merely a faint extended smudge when imaged with Keck 
in the R band for « 600 seconds. In either case, it was 
not possible to obtain spectra. While all of the previous 
analyses in this paper were conducted as if these fields 
did not exist, we actually must determine a method of 
including them in the sample in order for our survey to 
be considered complete. Assuming that the fields were 
truly empty, i.e. the optical source was fainter than our 
magnitude limit on Keck, and that pointing errors or some 
other experimental errors did not result in radio-source 
misidentification, we can make two extreme, simplifying 
assumptions. On the one hand, we can assume that no 
DLAs are present toward these optically faint quasars and 
calculate a lower limit on Qg by including some average 
redshift path length for each object, where the average 
pathlength is determined from the known quasars in our 
survey. On the other extreme we can assume that each EF 
is actually empty because of the presence of a high column 
density, dusty DLA. We can assume that each EF contains 
an average to high column-density DLA and estimate an 
upper limit on Vlg. 
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Fig. 6. — R band images of the empty or extended fields. Expo- 
sures times are given in Table 3. The circle radius is fa 5 arcsec. 
All image orientations are the standard North (up). East (left). 



Using the results of the combined sample, we determine 
the average redshift path per quasar to be < AX > = 
2.66. Assuming that each of the eight EFs would con- 
tribute this amount of redshift path gives a new total red- 
shift path searched of AX = 338.43. The lower limit on 
fig, assuming that none of the EFs contained a DLA, is 

cylower 



1.08 X 10~ . The upper limit, assuming that 
each of the EFs contains a DLA of average column den- 
sity in our survey, A'^hi = 7.86 x lO^^cm"^, results in 
an upper hmit of O^pp'^'" = 1.41 x lO'^. While both the 
lower and upper limit on fig, derived by including the 
EFs, are clearly within the error of the combined value 
of Qg — 1.15^Q'3g X lO^'^, it is notable that the upper limit 
is only w 22% larger than fig, i.e. even if each EF contains 
a dusty DLA of our average A^hi, the effect on fig would 
be relatively small. 

We can take the analysis one step further and allow the 
average value of A^hi to exceed 7.86 x 10^°cm~^. To de- 
termine the minimum average column density DLA that 
would affect our results we assume that each EF contains 
a quasar at our average redshift and a DLA of fixed col- 
umn density which we vary from the lower limit, Ajji = 



10 
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Fig. 7. — This figure demonstrates the potential iinpaet that 
the empty fields could have on Qg, the neutral gas mass density. 
We calculate the resultant value of Qg assuming that each of the 
eight empty fields contains a DLA of the stated A^hi- The diamond 
with Icr error bars is the value of Qg found from the combined 
radio sample, calculated by ignoring the empty fields. The asterisk 
represents the value of Qg if each of the empty fields contained a 
DLA of our average Ajji = 7.86 X 10^'' cm^^. For reference, the 
optically selected SDSS survey value is plotted as an unfilled circle. 



^^ cm ^, to the generally observed highest column 



1 X 10 cm "'. We plot the results in 



2 X 10 

densities of iVni 
Figure 7. 

As previously stated, if we assume that each EF contains 
a DLA of average column density in our survey, A^hi — 
7.86 X 10^°cm~^, we derive an ilg = 1.41 x 10~^, indicated 
in Figure 7 by the blue asterisk. Compare this with the 
red point and error bar, our radio-selected survey value 
of fig, derived by ignoring the EFs. For reference, the 
optically selected SDSS-DR3 survey, at ng = 0.82^°;°^ x 
10^"^, is plotted in green. From Figure 7 it is seen that if 
each EF contained a DLA of Ahi ~ 10^^'^ cm~^ or larger, 
a relatively large value occurring in only « 15% of our 
sample, the impact of the EFs would be large enough to 
increase the resultant value of fig by ss 50%. 

There are, however, several arguments for why these 
EFs are most likely not Zem > 2 quasars extinguished by 
very dusty, high column density DLAs. If we extrapolate 
the H I frequency distribution function, /hi(A^, X) result- 
ing from the radio-selected sample, as seen in Figure 3, we 
would expect not more than 1 DLA with Ahi > 10^^ cm~^. 
If two or more high column density systems existed, the 
resulting /hi(A^, X) would be unphysical assuming galax- 
ies have declining surface density profiles. In this case, 
we would require a bimodal population consisting of high 
column density, high dust-to-gas ratio systems, such as 
molecular clouds, that would be missed in optical surveys. 

While we cannot rule out the existence of a bimodal 
population, we can determine exactly how many high col- 
umn density systems the current radio-selected distribu- 
tion function would predict. We plot the cumulative num- 
ber of DLAs above a certain minimum H I column density 
and extrapolate using our single power-law fit. From this 
plot, in Figure 8, it is apparent that we would expect only 



0.3 DLAs with Ahi > 10^^ cm"^. Additionally, in the case 
of the five fields containing faint extended emission, we can 
use scaling arguments to show that if we assume this emis- 
sion is actually the resolved quasar host galaxy, then the 
quasar would have to be a low redshift object and would 
not have been included in our Zem > 2.0 survey. Adopting 
the typical high-z quasar host galaxy scale length of ~ 12 
kpc (Kuhlbrodt et al. 2005), we can estimate the angular 
size at z = 2.0 to be w 1.4 arcsec. Careful inspection of 
the extended fields reveals extended blobs on the order of 
3 arcsec or larger, making these low z quasars that would 
not have been included in our survey. 

In an effort to determine the nature of the EFs, we are 
currently conducting an observing program on the Green 
Bank Telescope (GBT) to search for 21 cm absorption 
along the sightlines to these EFs. We will carry out a 
redshift path search from z = [0.5, 3.9] using the frequen- 
cies of « 300 MHz - 900 MHz. If a high column density, 
dusty system does exist along the line of sight and is block- 
ing out the quasar light, it is likely that we will detect it 
in absorption. 

5. ANALYSIS & DISCUSSION 

The best way of determining the significance of dust ob- 
scuration of quasars is to compare the results of magnitude 
limited and radio-selected quasar surveys for DLAs. The 
primary problem with this method has so far been the 
limited survey size of the radio-selected surveys and the 
resultant large error bars that preclude conclusive results. 
While the UCSD survey itself was slightly smaller than 
the previously published CORALS survey, combining the 
two surveys in effect doubles the size. However, the uncer- 
tainties of the combined sample, whose size is still more 
than an order of magnitude smaller than the current opti- 
cal samples, are so large that definitive statements remain 
elusive. 
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Fig. 8. — The cumulative distribution of the number of DLAs 
with a specific minimum A^hi- Ovcrplotted is the single power-law 
fit to the distribution function (dotted line). It is seen that only 0.3 
DLAs are expected with Ahi > 10^^ cm^ . 



To make a simple estimate of the number of radio se- 
lected DLAs necessary to conclusively answer the ques- 
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Fig. 9. — The results of bootstrapping error on samples of 
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tion of dust bias, we performed the following analysis. If 
we desire a result with errors of no larger than « 10%, 
we can perform a bootstrap error evaluation on a ran- 
dom sample of DLAs each time increasing the number of 
DLAs to determine how many DLAs are necessary to give 
the desired precision. In Figure 9 we plot the results of 
our bootstrap error estimation, normalized to one and off- 
set for each different sample type. The estimation was 
performed on random Nm samples with minimum A^hi — 
lO^*^"^ cm^^, and three different maximum column densi- 
ties, A^Hi — lO^^'^^cm"^, to match the upper limit of the 
combined radio sample, A'hi = lO^^'^^cm"^, and A^hi = 



10 cm^ . It can be seen that in a sample with a max- 
imum column density similar to the combined sample, the 
desired « 10% error is acheived with a sample of w 100 
DLAs. As the maximum column density is increased, the 
error bars increase as well. Note that the bootstrap errors 
of our actual sample (number of DLAs = 26, error « 30%) 
are slightly bigger than those of the randomly generated 
samples due to the fact that our sample contains only a 
few high column density systems. 

To definitively answer the question of dust bias we would 
ideally hope to at least approach the total redshift path 
searched by optical surveys in order to make valid com- 
parisons. In the combined radio-selected sample, the total 
cosmological redshift path searched was AX = 317.11. 
Compared with the total redshift path surveyed by the 
latest large optical survey, SDSS-DR3 with AX — 7333.2, 
we are still more than an order of magnitude smaller. 

The combined radio-selected central value of Jig is slightly 
higher than all of the optically-selected values, as plot- 
ted in Figure 5. However, when considering the la lower 
limits of the radio-selected values of ilg , no difference be- 
tween the magnitude-limited sample and the radio-selected 
samples can be ascertained. If we ignore the possibility 
of evolution in ^g we can compare the entire SDSS-DR3 



optically-selected survey over the complete redshift range 
with the radio-selected value. Excluding the Peroux point, 
the SDSS-DR3 value taken over one redshift bin, z = [2.2, 



5.5], gives ilg = 0.82_o;og xlO"^. Comparing this value 
with the la lower limit of the combined radio-selected 
value of n'-g'"" = 0.77 x lO'^ over the range z = [1.51, 
4.4], we see excellent agreement. This agreement in J7g 
between the radio and optically-selected surveys for DLAs 
is the best evidence for our conclusion that dust bias does 
not have a major effect on the results of optically-selected 
surveys. 
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